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INT:sCOUCTION 


The stress distribution in an incomplete tore loaded as shown in 
Pig. 1 ie of particular interest since it very closely approximates that 
in heavy close-coiled halical springs under axial teneion or compression, 
Necessarily the spring helix angle must be small, which is the case in a 
close-ceiled spring. By a heevy spring is meant one whese retio of mean 
diameter to cross-sectional diameter is such a value that the curvature 
ef the section must be considered. | 

It should be noted that the 
stress distribution arising from 


the losding in Fig 1. is not pure 


torsion in the usual sense, but Cz 

is a combination of torsion and 

direct shear. The probiem there- : 
fore resolves itself into one of Pig. ls 


finding a single stress function which defines the true stress distribution 
in a cross-sectien of the circular ring sector. 

Several solutions to the problem are in the literature, all of which 
by various means selve the differertial ecquetion arising from the conditions 
of compatibility. The first, by Michell (1) in 1899, used polynomial 
stress funetions and obteined solutions for approximately circular cross- 
sections. Goéhner (2) used successive approximations te appreach an exact 
solution. Shepherd (3) used a methed similar to both Géhner and Michell 
by finding a sequence of functions for apprcximately circular cross-sections 


and combining them linearly in such a manner thet the sum was a soluticn. 
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Wahl (4) obtained a selution using curved bar theory and assuming a dis~ 
placement of the center of rotation. Southwell (5) presented a formal sol- 
ution for an arbitrary cross-section with a view towards a "relaxation" 
approach. Frieberger (6) has presented an exact sclution for a circular 
crose-section by finding a stress function analogous to the ordinary torsion 
function and soiving the problem in toroidal harmonics. 

In this paper an approximate solution is obtained using the principle 
of least work, A stress function is found satisfying the equations of ecuil- 
ibrium and the boundary conditions and whese corresponding stresses make the 
strain energy a minimum. The solution of the differential equation of conpst—- 
ibility has therefore been replaced by the problem of minimizing the strain 
energy. In the energy method, the condition of minimum strain energy is 
equivalent to satisfying compatibility net in a point by point sense, but 
“en the average" throughout the body. 

The purpese of this investigation has been to answer two questions in 
the euthor's mind. Namely, in view of the fact that nowhere wae the author 
able toe find the energy method used in the literature: 

(1) Can the problem be solved by this method, and how do the results 
compare with those of other solutions? 

(2) Does the problem particularly lend itself to solution by 
energy mathods? 

It was found that the problem is not adaptable to an exact solution by 
energy methods, but by making some approximations, excellent results ars 


obtained that agree very closely with Frieberger's exact solutien. 
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PORMULATICN OF TH PROBLEM 





We will consider a sector of a circular ring with mean radius of curva- 
ture R and crosa=-sectional radius a. A lead P is applied to one terminal 
eress~section as shown in Fig. 2, the other remaining fixed. Cylindrical 
coordinates are used, where the 3 axis coincides with the toroidal axis, 
and the axis of the ring sector liee in the ré plane. 





Fig. 2. 


& increases positively as shown in the figure and r increases outward 
from the toroidal axis. Later in the solution the coordinates will be trans- 
formed, but for the present purpose of establishing a stress function satis- 
fying the equations of equilibrium, cylindrical ceordinates are most. convenient. 
Assuming sero body forces, from THEORY OF SLASTICITY, Timeshenke & Goedier, 
Equations (170) the differential equations of equilibrium are 
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Using the same assumptions made by Gthner in this case, namely that the 


enly non-vanishing stresses are Tro 4 Tze 2nd that the stress distribution in 


any cross-section is independent of @ these reduce to 


Tre 3Tze ai e 
Cr i ial 


This may also be written 


} 
[3 (Te) + 2 Te)| ake) 
A stress function & satisfying the above is 
GES a8 : Cotes GR’ S = 





Mie = rv Tia 
Where G is a constant (actually the modulus of rigidity). 
Therefore the stresses may be expressed as 
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At this point it is convenient to transform the cylindrical ceerdinates 
rs @ into toroidal coordinates e 5 © (refer to Fig. 3) 
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re R- pcos? 
If 4 is a function of r and g, where 


& e e 
ee e sin rem Fig. 3 
Then 
a4, 29 ar, a4 dz 
3e St 3p * 32 3p Sor oe 2E <psinp 
a? , ah ar , 34 37 Ge . 
ay or oY 32 Sy 22 =sinY © eet 
Substituting 
a¢ oF 


3° 2 aye (- cos?) + 22.(sin¥) 
a 
= = 3€ (psind) + ae Cocos #) 


Selving for ae and ae 
or oz 
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& _ (cosW) 2b a¥) 3a¢ 
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In a plane cress-section determined by § a constant 


‘se Tbe * - Th, cos? + Tzg Sin? 
a. 


Tyo? Trp sinh + Tz cos? 


Using Equations (2), (3) and (4) the follewing result is obtained. 
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The latter expressions relate the stress function and the etresses in 
the new system cf coordinates. 

It follows that since the shear stress Tee is nermal to the boundary, 
it must vanish everywhere on the boundary. This is true because the surface 
ef the body is free from any external forces. Using this cendition with 
Equation (5), it is apparent that 2° =0 and ® must be constant en the 
boundary. 

The circular ring sector we are considering is a singly connected body, 
hence the constant may be chosen arbitrarily. Therefore tie boundary cendition 
is taken as $ = co) everywhere on the boundary. 

The only action en a cross-section is a force P directed aleng the 
toroidal axis. This may be resolved into a force and a couple as shown in 


Fig. ‘he 
R is equivalent te 


Fig. kh. 


It is now seen that the two conditions of static equilibrium to be 
satisfied are that the resultant stress on a cross-section preduce 2a force P 
directed along the zg axis end a moment about the center PR. Theee requirements 
may be written a ay, 
re Lf [ Cposine + Tre cost) p ded 
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The strain energy per unit angle 8 is 


a aw 


(7)———- | Jf a J] (% + Tye )(R- pcos?) o do dy 


The method of solutien will now be to take the stress function in the form 
p: = “te file » where >; are suitably selected functions of OP and 
p , each of which satisfies the boundary condition 4: -©O when P5™ , 
The coefficients % are constants which are evaluated from the minimum 


condition of strain energy. 
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FIRST APPROXIMATION 





For a first approximation we shall take a function $ » satisfying 
the boundary condition that it vanish everywhere on the doundary, in the 
form $e (phot) (me + = Cos + ) e The reasons for this 
particular choice are discussed in Appendix 4. Taking the partial deriva- 
tives of © with respect to the two veriables p and + 


ae = Zea, + a, (3% 07) cosy and a “i (pho) sy 
2° R mae R 


Substituting in Equations (5), the fellewing expressions are 
obtained for te and Tyo. 


(p-0') apes 
i+, a pS amt aot Tye: oe -| Zea silaphe coy 


1 
The appearance of the term fe pcos?) in the stress equations 


makes the integration required in (6) and (7) very complicated and the 
results largely unmanageable in the evaluation of the unknown coefficients 
in  , (See Appendix 3). This is particularly true when additional terns 
are used in p for a higher order of approximation, and in the evaluation 
of the strain energy where the stresses appear as squared terms. 


Since £ is — less than unity, we may write 


co (ele cost) as a(f)cos+ + 3(Z)eos+ ee 


Utilising this expansion, the exact stress expressions (5) may be 
approximated as follows 


pthc atl 22 .,| 
. Tyas 6 [C14 2 Res) 28 ena, = 3] 


This particular form of approximation accomplishes the desired result 
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0. 
of limiting the highest power to which the ratios R and £ appear in the 


stress equations. 


Since $.: (por) and co © lp eae Sa 
The partial derivatives are 

ae, a ae | _ (36a) Cos p 

2e P =) R 

a®, ro Seley sin 

av ; ze: 


Substituting, we arrive at the following approximate expressions for 


Tee = G tees te) 510 ¥| 
Aig * o[ 2p. $ phtes 3p sno og 


Substituting these values of The and Tye in the first of Equations (6), 


the stresses. 
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and integrating we obtain 


4 
Gra ry a PR 
P= ae 2 ee Sue) S Gtra3 


The same result is obtained from the second condition of Equations (6). 
It follows that %, is fixed by the requirements of static equilibrium 
and o& , may now be determined by the condition of minimum strain energy 
aU 


eG, oye 
that =a, 


From See 7 A 


ak 
Y= EL Cunt typ Bele pordptpsy 
Substituting the stresses from Equation (9) and integrating 
au. Ra [/1a Lot 
a4, G \( es + G %)s| 


Setting a = © and solving for ~, 
a 
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Using these results in Equations (9) we arrive at the expressions 
for the first approximation of the stress distribution in a cress~section 


of the incomplete tore 


——— PR ie P22) and) 
4 R 


pe 1a" 
(10) 
Tye * = | Ie casy| 
At the point of maximum stress where Pra and Y=zo the above 
reduce to 
Tpe *° 
(12) 


[re] Sel £8) 


It is interesting to note at this point that for this particular soluticn, 
one of the unknown coefficients in ’ is determined directly from the require- 
ments of static equilibrium, and the other directly from the minimum strain 


energy condition without constraint arising from static equilibriun, 
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SECOND APPROXIMATION 





A closer approximation to the true stress conditions will result 
if higher order terms of a suitable nature are used in the streas function. 
We shall now take $ as 


2 
d. (erate t OP cos + e “28 cos*+ + 1a 8 sinty + 4S} 


Reasons for this particular choice of functions are discussed in Appendix A. 
Again employing the binomial expansion of = 7 we write approximate 


(R-pcos 
expressions for Toe and Tye . 


i ${( + 2.-£c0s » +30 05") 220, re (( rL-Rcos¥) ad, + 284, ct 2954, + a] 
® sa Rcos #432, Cos npn, t (irafeosy) By, + Be hal 2 4, ‘ 2fu,( 


Where 
> a (oe o®) 0a6 e 2-0 ty + : ar( = 
® oP cosy 7 : tee Sin* 


This is an extension of the device used befcre to limit the highest 
power to which the ratios = and + appear in each term of the stress 
equations. Since = and " occur in a like manner in 5 $, and , 
these latter terms are grouped together and treated in similar fashion when 
introduced into the approximate expressions for the stresses. 

Taking the partial derivatives, substituting and rearranging the terms 
for convenient integration, the following approximate expressions for ‘yo 
and Ty, are obtained. 
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From the first of the static equilibrium conditions in (6) (that the 
resultant stress must produce a force FP in the Z direction) it is again 
found that O° _» 

‘ene * 

The second static ecuilibrium condition (that the resultant stress 

must produce a moment about the center equal to PR) gives the following 


result. 
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However, since 1 , PR 
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Since ©, 5%. ,&%3 and %, will ultimately all centain the factor 
& ot some simplification of the algebra will be efforded if we make the 
following substitutions 

ee za, (Gane where N= 1525354 


Finally the constraining function derived from the conditions of static 


equilibrium to be used in minimizing the strain energy is 
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The work invelved in obtaining the partial derivatives of the strain 
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energy with respect to the unknown coefficients aw, > ~, yy and “4 will 
be aimplified by differentiating under the integral sign. 


Therefore a un 


i Ef [ek + De = £ cos? )p de dp 


2e is not required since &, has already been evaluated from the 
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conditions of static equilibrium. Since Qu : (Constant) en ; 
n Ls) 


it is convenient here to take = - After substituting for Tp b 


and Tt wo from Equations (12) and performing the required integration, 
the following expressions are obtained. 
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To minimize the strain energy and evaluate the unknown coefficients 
6, J b PB: Bs and 6 q » the method of Lagrangian multipliers will be 
used with the constraining function (13) established by the requirements of 
static equilibrium. The constant is appearing in the partial derivatives 
of the strain energy will be incorporated in the multiplier. Letting dN be 
a Lagrangian multiplier, and { (8, 4 B a B, A 6.) = © the constraining function, 


we may write 
= 4. 
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where Ne 1525354 
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we arrive at the following set of equatioss, the solution ef which will 
evaluate 3, ? 8, 3 6, $ B, and determine the stress function. 
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Using these results in Equations (12), we may now write expressions 
representing a second approximation of the stress distribution in a cross- 
section of the incomplete tore, 


in in2?| anere 
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(16) 
Myer FSCS 25) FG ase + (4-WerG)alers IG-a9p- (5s) tap 
At the point of maximum stress, where pro and «=P=zo , the above 
reduce to 
ters 
(17) 
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As is apparent from the foregoing development, further approximations 
utilizing additional terms in the stress function will result in extremely 
long and tedious calculations, This in itself is a limitation of this methed. 
Therefore at this point, assuming the solution to be a rapidly converging one, 
we will stop and intreduce actual values cf the ratio of cross-sectional 


radius to the mean radius cf curvature of tne tore in order te compare results 


with other solutions. 
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RESULTS 


The distribution of shearing stress on a horizontal diameter is shown in 
Fig. 5 and the circumferential stress distribution in Fig. 6. In both cases a 
ratio cf R/a equal to 4 is used since this realizes the worst condition (i.e. 
greetest curvature for a given cross-secticn) of any practical significance. 


The quantity S appearing as the ordinate in both curves is a dimensionless 
Tre ta, 
i 

diameter and the periphery. A stress distribution representing pure torsion 


quantity and is equal to , Since Toe vanishes on beth a horizontal 
of a straight circular shaft is shown by a dotted lines in both figures. It is 
seen that the maximum stress actually existing in the tore is considerably 
greater than that derived from ordinary torsion theory. Both curves are in 
good agreement with similar ones derived from the exact solution by Frieberger. 
Points from Frieberger's curves appear as the small circles in Fig. 5 and 6. 

In comparing the results cf this solution with others, namely Gthner, Wahl, 
and Frieberger, the point of maximum stress will be used as a reference with 


different values of /a. Table 1 gives values cf K in the expression 





. LER 
Wi | iesus? To:  K] for the several sclutions. 
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Table 1. indicates that the energy method applied to this problem 
preduces results which compare favorabiy with other selutions. It also 
appeers that the solution cenverges rapidly, since only five terns were 


used in the stress function. 


17 















— 


~* 





SS . 
& 







—_ Ss ae lla " 
oe en ee 
Pees so Gees <m t= = 
(i+ as @ hom i « ~— a A 


a .. — o— ae = 





‘seat ht Om: __ a ~*~ o.. 
a (§ —a, Gut - @ 
’ eo a *_— - _ 











Fig. 5. 


Distibution of shearing stress on a horizgntal 
diameter for a/R=1/4. (420,17 ). S=(Tye) (ta )/P. 
Frieberger's points are indicated by small circles. 
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Circumferential stress distribution for 
a/R=1/4. S=(Tyg)(1a“)/P. Frieberger's points 
are indicated by small circles. 





In discussing any conclusicns from this investigation, it would be 
appropriate to recall the two questions that prompted it. 

(1) Can the problem be solved by this methed, and how de the results 
compare with these ef cther solutions? 

(2) Dees the problem particularly lend itself to solution by energy 
methods? 

First, the method will work and acceptable results are obtained with a 
relatively few terms in the stress function. This is in itself worthy of note, 
since it allows a very complex problem tc be attacked by the more elementary 
metheds of mathematics. 

However, in reference te the second question, there are limitations both 
inherent in the energy method end peculiar to this particular application, that 
strongly indicate the problem ia not especially adapted to a solutien by energy 
methods. 

The energy method, except in unusual circumstances, does not provide an 
exact solution. Consequently, in the absence of an exact solution, there is no 
real basis for judging the results. The fact also that the energy method requires 
minimizing an integral, which is done only with extreme difficulty with any 
number of terms in the stress function, 4s a limitation to its adaptability. 

In conclusion then, it may be said thet this solution has the value of 
arriving at very geod results using a relatively uncompliceted stress function 
of only five terns. 
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AFr RNDIX A 


According to the principle of least work which is used in this solution, 
an exact atress function would require selecting from all functions that satisfy 
the boundary condition those which minimize the strain energy. 

Since in general this procedure is too difficult, a limited number W of 
suitable functions was selected to determine an oppeoninnts stress function. 

In choosing functions of © and in P5 2 Xi Pz » the 


L*0 
first consideration was the boundary condition $ =O when ido 


This condition was satisfied by taking each &, to contain the factor Ce* E) 
A A 
Then = (Gee!) Z. Mise (e, p) 
Leo 
With rectangular coordinates (F.4q) in mind, where & = P cos'P and 
ie Psin? » the next logical step was to express = S. (& p 1) in a power 
series, The first six terms of such a series were considered, namely those 


involving 


se, | ae 


Since a horizontal diameter ( Neo ) on a plane cross-section (@ is a 
constant) is an axis of symmetry for the & surface, } must be even in a") and 
not. contain terms involving odd powers of i « In general r will appear to all 
powers since ( = =o ) is not axis of symmetry. Therefore the remaining terms 


expressed as functions of 9 and | are 
2 2 2 L 
ie eces i , Carer 


cS 
The term a in ?, while not consistent with this line of reasoning, 
appeared as a result of the binomial expansion used in approximating the stresses. 


It was extracted from G8hners solution where a like approximation was used. 
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Consequently the stress function was taken in the form 


& ‘ (prot) Bo a, (£)cos + + a, (2eosty +4, (&)siaty + 4,3 


Here P was replaced by £ so that %; would in all cases be 


PR 
G@na* : 


In the first approximation the first two terms were used and in the second 


the product cf a dimensionless number and the factor 


all five were introduced in eo i 
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APPENDIX B 


| 
The appearance cf the term (R-p cosh)” in the exact equations relating 


the stresses a the stress function gives rise to the occurrence of integrals 


ar 
of the type aff ing ai ae arsint dost dod ns [ [a sunt cos¥ dod? 
CR- pcos?) os CR- pcos +)3 
evaluating the strain energy and in consideration of the conditions of static 


equilibrium. 


Taking the simplest form of the first type, where m=zi , ne oO and 6° S 


we have _pdpdt 
{ {ees i 


Integrating first with respect te and setting Rec and a eal b 
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Multiplying by AY and is 


fs dpep- sav a | e ot ge ff AP 
C+bcosP Crb cos? (c+b cos ¥\~ 


ao 6p oe i § ei 
etbcosp = ChE IL C+ bcos} eb] c+ boos Y 








Cos! (bss cos 


1 w~ 
| where C” yb 


} 
¥ct-b* 








: ’ : j & ll - 
artiades anoidauye demas wt at “/y).59 5) mead odd 20 mn aT 
, hired ee eee $e 
ok W bab WF ea P arz “9 \ — pk gb 9Feo2v “azo nan ahaa ow , 


e( e606 - 5% 5 6 6 *(¥ee2G-AD 
elsade Yo eoeidliwao wis Ye coldetebleace at bre yy At08N8 akeise oid 


* > _ _ ~~ |) —_ * 


° 28 bw O20 , leet omedw , ogg Jerlt edd to —s ot: 
& 

Pha ~ my 

“4 z02Q - -A) S 

d=G- pe 232A guttgec bus Y 62 Soequai ddbw satlt gntie 


~~" % 
“Yveoodt>) 


riz ) 
“Cpeood +o) © a 


weosotd , Wieer-')d + (w cords 2)¥ 209 | Ib 
“(20rd +>) “(yeood+>) Yb 








“Weord +9) yreod to) d “(py 200 d+) 


gritergesat bos Pb ye uty fot oN 


Poza |= -( iz _ (veord 49) F 4 G4 : 


4 
° 


_ ¥Jo_ $3. Sof) “> Wd =a aie Gg = tb ¢ = 
~& 200 d 02 d+) “da. eood+ald ~ Peord+o 
oe) 3 ,f =e -. vb 
® 2ood +) “dp Peoad+ DOD} YD ~ Peoddyd 





= t 
WY zor D+ id : 2oD 


~ s 
_ — tea %a~"D 


Therefore 


Ae ae Gree % c ee b+ e cos? 
(c+ bcos?) cbr \Cr bcos (c- By C+b cos? 


Intreducing the limits © and 21 , this reduces to 





21e C=zR 
(ch pryv2 where ig = e 


Therefore at 
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The other more complicated forms where nfo and Lit © are integrable 
in finite terms by similar reduction methods, out it is apparent that the work 
becomes excessively involved. Also the results in the form just developed are 
not readily usable in evaluating the unknown coefficients in the stress function. 

In view of the foregoing, despite the fact that it was not actually necessary, 
it was expedient te approximate the stresses in such a man:.er that the integration 
was simplified and the results put in a usable forn. 

This device of approximating the stress equations compromised the requirement 
that the stresses satisfy the equations of equilibrium. However, it appears that, 
since the stresses do satisfy the conditions cf minimum strain energy and stetic 


equilibrium, and give satisfactory results, the compromise may be tolerated. 
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